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The NRQCD factorization approach for calculating inclusive production of heavy quarko- 
nium gives unambiguous predictions for the polarization of quarkonium states. The 
factorization formula for polarized states can be obtained by using the threshold ex- 
pansion method to calculate the short-distance coefficients and then using symmetries of 
NRQD to reduce the NRQCD matrix elements. A particularly dramatic prediction of the 
NRQCD factorization framework is that prompt psi's and psi-primes's produced at the 
Tevatron should be predominantly transversely polarized at large transverse momentum. 

1. Introduction 

The NRQCD factorization approach0 is a systematic framework for analyzing the 
inclusive cross sections and annihilation decay rates of heavy quarkonium. The 
cross section (or decay rate) is factored into short-distance coefficients that are 
computable using perturbation theory and long-distance NRQCD matrix elements. 
The matrix elements scale as a definite power of v, the typical relative velocity of the 
heavy quark in quarkonium. Thus the cross section can be organized into a double 
expansion in powers of v and powers of asirric), where tUc is the heavy quark mass. 

The NRQCD factorization framework makes definite, and in some cases rather 
dramatic, predictions for the dependence of the cross section on the polarization of 
the quarkonium. Below, 1 outline the NRQCD factorization approach as it applies 
to polarized quarkonium states and then summarize the applications that have been 
carried out thus far. 



2. NRQCD Factorization 

The general expression for the cross section for the inclusive production of a quarko- 
nium state H with four-momentum P is 

1 P 
y"da(12^ H(P) + X) = 

*Invited talk presented at the Quarkonium Physics Workshop, University of Illinois at Chicago, 
June 1996. 
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X ^ {2T:f5\h + k2-P-kx) \T,2^HiP)+xf , (1) 

X 

where Ti2^h{p)+x is the T-matrix element and the sum on the right side includes 
integration over the phase space of the additional particles. This cross section 
involves both "short distances" of order l/rric or smaller and "long distances" of 
order l/{mcv) or larger. The production of the cc pair involves short distances, 
because the parton processes that produce a cc pair always involve particles that 
are off their mass shells by amounts of order rric and can therefore propagate only 
over short distances. The binding of the c and c into the state H involves long 
distances, because gluons whose wavelengths are comparable to or larger than the 
size of the bound state, which is of order l/(mcv), play a large role in the binding. 

If the cross sectionin is sufficiently inclusive, short-distance and long-distance 
effects can be factoredlll. (If there are hadrons in the initial state, it is also necessary 
to restrict the four-momentum of the quarkonium to be significantly different from 
that of the initial hadrons in order to avoid contributions from diffractive scattering.) 
Like many other factorization "theorems" of perturbative QCD, the factorization 
of quarkonium cross sections has not been proven with complete rigor, but it is a 
plausible generalization of the factorization theorems for the Drell-Yan production of 
muon pairs and for heavy quark production. The factorization relies on cancellations 
between soft partons that are emitted by the cc pair and soft partons that are 
exchanged between the cc pair and other jet-like collections of coUinear partons. 
After taking into account these cancellations, the cross section can be factored 
into short-distance and long-distance parts. The short-distance parts involve the 
production of a cc pair with small relative momentum plus hard partons. The cc 
pair that emerges from the short-distance part is essentially pointlike on the scale 
of the quarkonium wavefunction. The long-distance parts involve the formation of 
H plus soft partons from the pointlike cc pair. 

The standard factorization methods of perturbative QCD produce an expression 
for the cross section that involves an integral over the relative momentum of the c 
and c that form the quarkonium. Long-distance and short-distance effects can be 
further untangled by expanding the short-distance factors in powers of the relative 
momentum q and absorbing the integration over q into the long-distance parts. The 
resulting long-distance factors can be expressed as matrix elements in an effective 
field theory called nonrelativistic QCD (NRQCD). The T-matrix elements in (|l|) 
that survive after soft-parton cancellations can be expressed in the form 

^12^ff(P)+X„+Xs ~ ^'il2~>cc{Pm)+XH{H + Xs\tp^lCnX{x ^0)\0) , (2) 
n 

where the sum includes all possible color and angular-momentum states of the 
cc pair. The factor Ti2^cc{p,n)+XH '^^^ interpreted as a T-matrix element for 
producing a cc pair in the state n plus the hard partons Xh ■ The second factor on 
the right is a matrix element in NRQCD between the vacuum state and a state that 
in the asymptotic future consists of the quarkonium H at rest plus the soft partons 
Xs- The local operator tp'^ICnX creates a pointlike cc pair in the state n. 
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After inserting the expression for the T-matrix elements into (|l|), we can 
rearrange the cross section into a form in which the short-distance and fong-distance 
contributions are factored: 

X 

~ X! ( y^X^l2^cc(P,m)+XH)*'^12^cc(P,n)+XH 1 
mn \ Xh / 

X (^{Olx^lClm + Xs){H + Xs\i^^lC,,x\0)^ ■ (3) 
The short-distance factor is 

Cmn{kl,k2,P) = ^(^2— cc(Km)+XH)*^2->cc(P,n)+Xff • (4) 
Xh 

The fong-distance factor is 

(OD - {0\x^K.l^VH^^1CnXm. (5) 

where Vh projects onto states that in the asymptotic future contain the quarkonium 
state H plus soft partons: 

Vh = Y.\H + Xs){H + Xs\. (6) 

Xs 

Inserting ^ into (|l]), we obtain the NRQCD factorization formula for the inclusive 
cross section: 

Y.daiU^H[P)+X) - (^^)3^E^™"(fci,fc2,P)(0-„). (7) 

X mn 

We can write a similar factorization formula for the inclusive decay rate of the 
quarkonium H via the annihilation of the cc pair. The general formula for the decay 
rate is 

Y,dT{H-.X) = ^ E {'2^t5\P-kx) \Th(p)^x\\ (8) 

X " X 

where P = [Mh, 0). This can be expressed in the factored form 

dm ^X) ^ ^rnn {H\0^n\H) , (9) 

X ^ mn 

where the Cmn's are short-distance coefficients. The long-distance factors are ex- 
pectation values in the quarkonium state of local four-fcrmion operators of the form 

0,nn = ^p^K.l^XX^'K^nlp- 
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3. Short-distance Coefficients 

Since the coefficients Cmn in (0) and (||) involve only short distances of order 1 / rric 
or larger, they can be expressed as perturbation series in as{mc). These coefficients 
are known at tree level for many processes, and in a few cases they are known at 
the one-loop level. Most of these coefficients have been obtained by calculating 
the perturbative cross section for producing a cc pair in a state with a prescribed 
nonrelativistic wavefunction. The inclusive cross section is sensitive only to the 
behavior of the wavefunction near the origin. The analogs of the NRQCD matrix 
elements (|^) for the cc state can also be calculated using perturbation theory in 
terms of the prescribed wavefunction. Knowing the cross section and the matrix 
elements, we can read off the short-distance coefficients. Unfortunately, this method 
is not sufficiently general to determine all the short-distance coefficients, especially 
for polarized quarkonium states. 

The threshold expansion method, developed recently by Braaten and Chera, is 
a general prescription for calculating the short-distance coefficients. It is based 
directly on the NRQCD factorization approach, and can be readily applied to po- 
larized quarkonium. 

3.1. Threshold expansion method 

The threshold expansion method relies on the fact that the short-distance coef- 
ficients in are insensitive to the long-distance effects that bind the cc pairs 
into the quarkonium state H. Thus the factorization formula will hold with the 
same short-distance coefficients if we replace H by asymptotic perturbative states 
cc = cc(q, ^, 77) that consist of a c and a c with relative momentum q and spin/color 
state that is represented by the Pauli spinors ^ and rj. To completely determine the 
short-distance coefficients, we need to use different states cc and cc' in the T-matrix 
element and in its complex conjugate. The resulting matching prescription is 



J2 i2Tr)^S\h +k2-P-kx) (ri2^ee'(P)+x)*ri 



12^cc(P)+X 



X 



pQCD 



« Vc™„(fci,fe,P)(o|x^/c„i/'Pec',ccV'^/c„x|o) , (10) 

^ pNRQCD 
mn 

where the projection operator in the matrix element is 

Vcc',cc = J2\cc' + Xs){cc + Xs\. (11) 

Xs 

The left side of ( p^ is to be calculated using perturbative QCD, and then expanded 
in powers of the relative momenta q and q'. The matrix elements on the right side 
are to be calculated using perturbative NRQCD, and then expanded in powers of 
q and q'. The coefficients Cmn are then determined by matching these expansions 
order by order in as- 



3.2. Example 
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We illustrate the threshold expansion method by carrying out one of the simplest 
matching calculations. It gives the short-distance coefBcient corresponding to the 
parton process qq cc. The T-matrix element for this process is 



Ti2^ce = vik2h^.T''u{k^) u(p)7'^T%(p). (12) 

Making a nonrelativistic expansion of the spinors of the c and c, this reduces to 

Ti2^cc - 2^ vik2h^T-u{h) L^^ ^^a'T'^rj, (13) 

where L^j are elements of the boost matrix that transforms from the rest frame of 
the cc pair to the frame in which it has total four-momentum P. Multiplying by the 
complex conjugate of 7i2~,cc' and averaging over initial spins and colors, we obtain 

/I 2 2 

where z is a unit vector in the direction of the momenta of the colliding c and c. 
The spinor factor can be expressed in terms of an NRQCD matrix element: 

(xV^T''V|cc')(cc|V'V*r"x) = W V^T'^C'^V'T'^r/ . (15) 

Using the matching prescription (p^), the short-distance coefficient of the matrix 
element (xVJT^VT'ce'.ccV'^cr'T^x) is 

a, = {27:)^d\ki + k2-P)^ [6^' - zH^] . (16) 

Inserting the short-distance coefficient (^^ into the factorization formula (0) 
and integrating over the phase space of the quarkonium, we get an expression for 
the inclusive cross section: 



X 

3 2 

= - ^rnl) ^ [S^^ - z^z^ (xV-'T> Vh V'V'T'^x) • (17) 

This is a term in the factorization formula for the inclusive cross section of any 
quarkonium state H , whether polarized or unpolarized. There are additional terms 
of order aj. from the process gg —> cc. All other parton processes give terms with 
short-distance coefficients of order or higher. 

4. NRQCD Matrix Elements 

The long-distance factors in the NRQCD factorization formulas are expressed as 
matrix elements of local four-fermion operators in NRQCD. Since long-distance 
effects in QCD are inherently nonperturbative, the NRQCD matrix elements can 
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only be calculated using nonperturbative methods like lattice gauge theory. There 
are effective lattice prescriptions for calculating the matrix elements {H\Omn\H) 
that appear in quarkonium decay ratesU. Unfortunately, these methods cannot be 
used to calculate directly the production matrix elements (0^„). The problem lies 
in implementing on the lattice the projection defined by (|^). In the absence of 
nonperturbative calculations, the only alternative is to treat the NRQCD matrix 
elements as phenomenological parameters to be determined by experiment. 

4.1. Model-independent framework 

The factorization formula provides a model-independent framework for analyz- 
ing quarkonium production. In any reasonable model for the production of quarko- 
nium through short-distance parton processes, the inclusive cross section can be 
expressed in the factored form (0). The model can therefore be reduced to a set 
of assumptions about the NRQCD matrix elements. Until recent years, most cal- 
culations of quarkonium production were carried out using either the color-singlet 
model or the color- evaporation mode@. In the color-singlet model, the quarkonium 
is assumed to be simply a color-singlet cc pair in an appropriate angular-momentum 
state. Only one NRQCD matrix element is assumed to be important, and it can be 
expressed in terms of the cc wavefunction, or one of its derivatives, evaluated at the 
origin. In the color-evaporation model, the color and angular-momentum quantum 
numbers of the quarkonium are simply ignored. The NRQCD matrix elements are 
assumed to be calculable in perturbation theory up to an overall normalization con- 
stant that depends on the state H . This model implies that the matrix elements 
scale like v^^^ , where w is a small parameter and D is the number of covariant 
derivatives D in the operator O^n. 

NRQCD predicts a much more intricate hierarchy among the matrix elements. 
The matrix element (C^„) defined in (|) scales like y'i+D+E+2M ^ where D is the 
number of covariant derivatives D that appear in the operator and E and M are 
the number of chromoelectric and chromomagnetic transitions that are required for 
cc pairs in the states created by ip^JCmX ^-nd ip'^JCnX to reach the dominant Fock 
state of H. These velocity-scaling rules determine the approximate magnitudes of 
NRQCD matrix elements. By keeping only those matrix elements that scale with 
the fewest powers of v, we can reduce their number sufficiently that a phenomeno- 
logical approach becomes tractable. One should keep in mind, however, that the 
importance of a term in the cross section is determined not only by the magnitude 
of the matrix element but also by the magnitude of its short-distance coefficient. 

4.2. Reducing the matrix elements 

The matrix elements can be further simplified by using symmetries of NRQCD. To 
illustrate the simplifications, we will use matrix elements of the J/V'(A), where the 
polarization state is specified by the helicity A. 



4.2.1. Rotational symmetry 
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Rotational symmetry is an exact symmetry of NRQCD. It implies, for example, 

that the matrix element for H = tpi^) in (O) must be a linear combination of the 

J. J. ' — 

tensors Sij , U\j , and U\j U-^ , where Ui\ is the unitary matrix that transforms 

vectors from the spherical basis to the Cartesian basis. If that matrix element is 

summed over the polarizations of the ip, the only possible tensor is 6ij. The matrix 

element must therefore satisfy 



4.2.2. Heavy-quark spin symmetry 

Heavy-quark spin symmetry is an approximate symmetry of NRQCD that holds up 
to corrections of order . The symmetry follows from the fact that the spin of the 
heavy quark is conserved at leading order in in NRQCD. It implies, for example, 
that the matrix element for H = ■0(A) in ( pT| ) must be proportional to U\jUl^: 

Summing over helicities, we recover (p^. 

4.2.3. Vacuum- saturation approximation 

The vacuum-saturation approximation can only be applied to specific color-singlet 
matrix elements. If the matrix element is expressed in the form (|^), the operators 
V'^^mX s-nd Tp'^lCnX must create pointlike cc pairs in the dominant Fock state of 
the quarkonium. In the vacuum-saturation approximation, the projection operator 
Ph defined in is replaced by \H){H\, which corresponds to keeping only the 
vacuum term in the sum over soft states Xs- This is a controlled approximation 
in NRQCD, holding up to corrections that are of order f^. The vacuum saturation 
approximation can be illustrated by the following matrix element for J /ijj produc- 
tion: 

{0\x^a'^V^(x)^^a^x\0) « (0|xVVI^(A)) (V'(A)|VV^x|0) . (20) 
The vacuum-saturation approximation can also be used for decay matrix elements: 

(^IxVv- • ^VxlV') « ^|(V'(A)|vVx|0)|'. (21) 

A 

The vacuum-to-quarkonium matrix element ('0|'0t(Tx|O) that appears on the right 
sides of ( pO| ) and (|2^) is proportional to the wavefunction at the origin. This matrix 
element can be easily calculated using lattice simulations of NRQCD. Thus the 
vacuum-saturation approximation provides a way to calculate certain production 
matrix elements on the lattice. 



5. Polarization Predictions 
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The NRQCD factorization formulas apply equally well to any quarkonium state iJ, 
including a polarized state. Since the short-distance coefficients are independent 
of iJ, the dependence on H enters only through the NRQCD matrix elements. In 
particular, the dependence on the polarization comes only from the matrix elements. 
In many cases, it is completely determined by the symmetries of NRQCD. 

The color evaporation model predicts that quarkonium states are always pro- 
duced unpolarized. Both the color-singlet model and the NRQCD factorization 
approach give nontrivial predictions for polarized quarkonium. There have been 
many calculations of polarization effects in the color-singlet model. In most cases, 
the resulting predictions are not particularly dramatic. Below, we discuss several 
examples for which the calculations have have been extended to include color-octet 
mechanisms predicted by the NRQCD factorization formalism. In one case, we find 
a very dramatic prediction. 

5.1. Spin alignment at the Tevatron 

The NRQCD factorization framework has led to a dramatic change in our under- 
standing ^f the production of charmonium at lar^^e transverse momentum in pp 
collisionsH. As pointed out by Braaten and YuanO in 1993, the cross section for 
pp ^ + X a,t sufficiently large transverse momentum pT is dominated by gluon 
fragmentation. It can be factored into the cross section for producing a gluon with 
large transverse momentum and a fragmentation function: 



The fragmentation function Dg^^{z) gives the probability that the jet initiated by 
the gluon includes a if) carrying a fraction z of the gluon momentum. Using the 
NRQCD factorization approach, the fragmentation function can be expressed in the 
form 



The matrix element that is leading order in v is \ {4'\4>^ crx\Q)\'^ , which scales like v'^ . 
It has a short-distance coefficient of order a^. Using this term in the fragmentation 
function (p3|), the cross section predicted by (^) is about a factor of SOielow recent 
data on prompt ip production at the Tevatron from the CDF detectorQ. 

In 1995, Braaten and Fleming! suggested that the gluon fragmentation function 
for the ip might actually be dominated by a term that represents a color-octet 
production mechanism. The matrix element is {x^a'^T'^tjj ^^a''T°'x)i which is 
of order v"^ and measures the probability of producing a ip from a pointlike cc pair in 
a color-singlet ■^S'l state. The reason this matrix element might be important is that 
its short-distance coefficient is of order as- The enhancement from the two fewer 
powers of as can overcome the suppression by u**. The leading-order expression for 





(23) 



mn 
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this term in the fragmentation function is 

D,^^ (z) ^ ^S{l-z) (x V^T'^V n V' V^T'^x) • (24) 

The pT-dependence predicted by this mechanism is in agreement with the CDF data. 
The normahzation depends on the unknown matrix element in (p^. The value of 
the matrix element required to fit the CDF data is consistent with suppression by a 
factor of relative to the corresponding color-singlet matrix element | (Tx|0) 

Cho and Wised pointed out in 1995 that this production mechanism has dra- 
matic implications for the polarization of the ip. At leading order in Og, the tp's 
produced by gluon fragmentation will be 100% transversfily polarized. The radia- 
tive corrections were examined by Beneke and Rothsteinll9, who concluded that the 
spin alignment at large pr will remain greater than 90%. The dominant corrections 
to the spin alignment at the values of px measured a,t the Tevatron come from 
nonfragmentation contributions to the cross sectionO. At large px, they fall like 
1/Px relative to gluon fragmentation. Unfortunately, the calculations required to 
obtain the prediction for the spin alignment as a function of px have not yet been 
carried out. 

The predictions for the spin alignment of the ip' are identical to those for the ijj. 
The spin alignment for the ■0' is easier to measure, because, in the case of the ip, 
one has to take into account the effects of the radiative decays Xcj '07- It should 
be possible to make at least a crude measurement of the spin alignment of the tp' 
from existing CDF data. 

5.2. Spin alignment in decay 

The color-singlet model predicts that the cross sections for producing charmonium 
in decay are too smaJito be observed at LEP. As pointed out by Cheung, Keung, 
and Yuan and by Choll3, the dominant contribution comes instead from a color- 
octet production mechanism involving the same matrix element that appears in the 
fragmentation function (|2j). Using the value of this matrix element obtained by 
fitting the CDF data, they found that the rate is almost an order of magnitude larger 
than predicted by the color-singlet model and thus large enough to be observed. The 
spin alignment of the %jj in Z'^ decay was also calculated by Cheung, Keung, and 
Yuan. Unfortunately, the alignment is predicted to be small and is completely 
unobservable in the present data sample from LEP. 

5.3. Spin alignment in B decay 

In the production of from B decay, color-singlet contributions are suppressed by 
a near cancellation between Wilson coefRcients in the effective weak hamiltonian. 
Color-octet production mechariiams are therefore important, in spite of the d"* sup- 
pression of the matrix elementscJ. The spin alignment of the i/j's that are produced 
in the decay B ^ tp + X were recently calculated by Fleming et al.lU. It depends 
sensitively on the values of 3 independent color-octet matrix elements, including 
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the one that appears in (p4|). A measurement of the spin ahgnment would therefore 
place strong constraints on these matrix elements. 

6. Conclusions 

The factorization formulas (|^) and provide a model-independent framework 
for analyzing heavy quarkonium production and annihilation rates. All the short- 
distance factors can be calculated systematically using the threshold expansion 
method. The long-distance factors are defined in terms of NRQCD matrix elements. 
The decay matrix elements can be computed using lattice simulations of NRQCD, 
but most of the production matrix elements must be treated as phenomenological 
parameters. The relative magnitudes of the matrix elements are predicted by the 
velocity-scaling rules of NRQCD. These magnitudes have a pattern that is very 
different from that assumed in the color-singlet model or in the color-evaporation 
model. 

The NRQCD factorization approach gives unambiguous predictions for the po- 
larization of heavy quarkonium states. These predictions are inescapable conse- 
quences of this framework. The polarization predictions from NRQCD factoriza- 
tion can be dramatic. An example is the spin alignment of prompt ip and tp' at 
the Tevatron, which is predicted to be greater than 90% at the largest values of pt- 
An experimental measurement of this spin alignment would be a crucial test of the 
NRQCD factorization framework. 

Acknowledgements 

This work was supported in part by the United States Department of Energy, Di- 
vision of High Energy Physics, under grant DE-FG02-91-ER40690. 

References 

1. G.T. Bodwin, E. Braaten, and G.P. Lepage, Phys. Rev 1)51, 1125 (199 5). 

2. E. Braaten and Y.-Q. Chen, OHSTPY-HEP-T-96-010 ( tiep-ph/9604237|). 



3. G.T. Bodwin, D.K. Sinclair, and S. Kim, ANL-HE P-PR-96-28 ( |hep-lat/9605023| ) 



4. G.A. Shuler, CERN-TH. 7170/94 ( tiep-ph/9403387| ), to appear in Ph ysics Reports. 



5. E. Braaten, S. Fleming, and T.C. Yuan, OHSTPY-HEP-T-96-001 (|hep-ph/960237^ ) 



to appear in Annual Reviews of Nuclear and Particle Science. 

6. E. Braaten and T.C. Yuan, Phys. Rev. Lett. 71, 1673 (1993). 

7. A. Sanson! et al. (CDF Collaboration), FERMILAB-CONF-95/263-E, to appear in the 
Proceedings of the Sixth International Symposium on Heavy Flavor Physics. 

8. E. Braaten and S. Fleming, Phys. Rev. Lett. 74, 3327 (1995). 

9. P. Cho and M. Wise, Phys. Lett. B346, 129 (1995). 

10. M. Beneke and I.Z. Rothstein, Phys. Lett. B372, 157 (1996). 

11. P. Cho and A.K. Leibovich, Phys. Rev. D53, 6203 (1996). 

12. K. Cheung, W.-Y. Keung, and T.C. Yuan, Phys. Rev. Lett. 76, 877 (1996); P. Cho, 
Phys. Lett. B368, 171 (1996). 

13. P. Ko, J. Lee, and H.S. Song, Phys. Rev. D53, 1409 (1996). 



14. S. Fleming et al., MADPH-96-953 (hep-ph/9608413) 



